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ABSTRACT 



Context. The frequent presence of weak magnetic fields on the surface of spotted late-B stars with HgMn peculiarity 
f\f ■ in binary systems has been controversial during the two last decades. Recent studies of magnetic fields in these stars 

using the least-squares deconvolution (LSD) technique have failed to detect magnetic fields, indicating an upper limit 
on the longitudinal field between 8 and 15 G. In these LSD studies, assumptions were made that all spectral lines are 
C"| identical in shape and can be described by a scaled mean profile. 

O [ Aims. We re-analyse the available spectropolarimetric material by applying the moment technique on spectral lines of 

inhomogeneously distributed elements separately. Furthermore, we present new determinations of the mean longitudinal 
magnetic field for the HgMn star HD 65949 and the hotter analog of HgMn stars, the PGa star HD 19400, using FORS 2 
installed at the VLT. We also give new measurements of the eclipsing system ARAur with a primary star of HgMn 
peculiarity which were obtained with the SOFIN spectropolarimeter installed at the Nordic Optical Telescope. 
Methods. We downloaded from the European Southern Observatory (ESO) archive the publically available HARPS 
spectra for eight HgMn stars and one normal and one superficially normal B-type star obtained in 2010. Out of this 
sample, three HgMn stars belong to spectroscopic double-lined systems. The application of the moment technique to 
the HARPS and SOFIN spectra allowed us to study the presence of the longitudinal magnetic field, the crossover effect, 
and quadratic magnetic fields. Results for the HgMn star HD 65949 and the PGa star HD 19400 are based on a linear 
I | regression analysis of low-resolution spectra obtained with FORS 2 in spectropolarimetric mode. 

. Results. Our measurements of the magnetic field with the moment technique using spectral lines of several elements 

■ separately reveal the presence of a weak longitudinal magnetic field, a quadratic magnetic field, and the crossover effect 
on the surface of several HgMn stars as well as normal and superficially normal B-type stars. Furthermore, our analysis 
suggests the existence of intriguing correlations between the strength of the magnetic field, abundance anomalies, and 
binary properties. The results are discussed in the context of possible mechanisms responsible for the development of 

■ the element patches and complex magnetic fields on the surface of late B-type stars. 

Key words, stars: chemically peculiar — stars: abundances — stars: atmospheres — stars: binaries: spectroscopic stars: 
individual (HD 11753, 41 Eri, 66Eri, ARAur, ftCnc) — stars: magnetic field — stars: variables: general 

X: 

H . 

C« ■ 1. Introduction troscopic binaries (SB) with late B-type primaries (spectral 

types B7-B9) with the goal of understanding why the vast 

The origin of abundance anomalies observed in late B-type ma j ority of these stars exhibit in their atmospheres cer- 

stars with HgMn peculiarity is still poorly understood. Over tain chemical abundance anomalies, i.e. large excesses of 

the last few years, we have performed extensive spectro- p Un> G&) Br SVj y , Zr, Rh, Pd, Xe, Pr, Yb, W, Re, 

scopic studies of both single late B-type stars and of spec- Q ^ pt; All) and Hgj and un derabundances of He, Al, Zn, 

Ni, and Co (e.g. Castelli & Hubrig 2004a). Strong isotopic 

* Based on observations obtained at the European Southern anomalies were detected for the chemical elements Ca, Pt, 

Observatory (ESO programmes 076.D-0169(A), 076.D-0172(A), and Hg with patterns changing from one star to the next 

^ 4 -M338(A), 085 £-0296(A), 085.D-0296(B), 087 D-0049(A), (Hubr f et al rfgggg £ Hubri 2004b; Cowl et 

088.D-0284(A , SOFIN observations at the 2.56 m Nordic \ nwrm r^u T 7 ^ n 4.1, — Z — t ■ 

„ ,. , „ , v t ni j u i , . j . , i al. 12010). Observationally, these stars are characterised by 

Optical telescope on La Palma, and observations obtained with . — ' . . . .,. /, . .,.„„, _i ., , , 

the CORALIE Echelle Spectrograph on the 1.2 m Euler Swiss lactational velocities ({v smi) < 29 kms , Abt et al. 

telescope on La Silla Chile I1972|) . The number of these chemically peculiar stars, usu- 
** F.R.S.-FNRS Postdoctoral Researcher, Belgium 



1 



S. Hubrig et al.: Magnetic fields of HgMn stars 



ally called HgMn stars, decreases with increasing rotational 
velocity (Wolff & Wolff IT974)) . 

More than two thirds of the HgMn stars are known to 
belong to spectroscopic binaries (Hubrig & Mathvs f!995|l . 
with a preference of orbital periods ranging from 3 to 20 
days. It is striking that the inspection of SB systems with a 
late B-type primary in the 9 th Catalogue of Spectroscopic 
Binary Orbits (Pourbaix et al. I2009[) indicates a strong 
correlation between the HgMn peculiarity and member- 
ship in a binary system among bright, well-studied SB sys- 
tems with late B-type, slowly rotating (v sini < 70 kms -1 ) 
primaries with an apparent magnitude of up to V ~ 7 
and orbital periods between 3 and 20 d. With the excep- 
tion of HR7241, all 21 systems have a primary with a 
HgMn peculiarity. Based on this fact, it is very likely that 
the majority of slowly rotating late B-type stars formed 
in binary systems with certain orbital parameters become 
HgMn stars. Consequently, careful studies of these pecu- 
liar stars are important for the general understanding of 
B-type star formation in binary systems. According to Abt 
& Snowden (1973), close-binary formation is inhibited when 
strong magnetic fields are present. They also suggest that 
close binaries dissipate the magnetic field. 

The aspect of inhomogeneous distribution of some 
chemical elements over the surface of HgMn stars was first 
discussed by Hubrig & Mathys (|1995p . From a survey of 
HgMn stars in close SBs, it was suggested that some chem- 
ical elements might be inhomogeneously distributed on the 
surface, with, in particular, preferential concentration of 
Hg along the equator. The first definitively identified spec- 
trum variability was reported for the binary HgMn star 
a And by Wahlgren et al. ([2001]) and Adelman et al. ([20021 . 
who showed that the spectral variations of the Hgn line 
at A3984 discovered in high-dispersion spectra are not due 
to the orbital motion of the companion, but produced by 
the combination of the 2.8 d period of rotation of the pri- 
mary and a non-uniform surface distribution of mercury, 
which is concentrated in the equatorial region. Their re- 
sults are in good correspondence with those of Hubrig & 
Mathys (fTM5| . The variability of the Hgn line at A3984 
was interpreted with a Doppler imaging code, revealing 
high-contrast mercury spots located along the rotational 
equator. Using a Doppler imaging reconstruction of spec- 
troscopic time series obtained over seven consecutive years, 
Kochukhov et al. (|2007p suggested the presence of a secular 
evolution of the mercury distribution. On the other hand, 
the work of other authors has proved that not only mercury 
abundance appears distributed in patches over the stellar 
surface. Almost all other elements, most typical Ti, Cr, Fe, 
Mn, Sr, Y, and Pt, are concentrated in spots of diverse size, 
and different elements exhibit different abundance distri- 
butions across the stellar surface (e.g. Hubrig et al. [2t)06a 
Briquet et al. 120101 Makaganiuk et al. I2011a|l . Moreover, 
an evolution of the element abundance spots at different 
time scales was discovered in two additional HgMn stars. 
Briquet et al. (|2010[) reported the presence of dynamical 
spot evolution over a couple of weeks for the SB1 system 
HD 11753, while Hubrig et al. (2010) detected a secular ele- 
ment evolution in the double- lined eclipsing binary AR Aur. 
Importantly, recent results (e.g. Nunez et al. 120111 Hubrig 
et al. I2011[) show that line-profile variability of various el- 
ements caused by a non-uniform abundance distribution is 
a general characteristic of HgMn stars, rather than an ex- 
ception. 



Typically, inhomogeneous chemical abundance distribu- 
tions are observed only on the surface of magnetic chem- 
ically peculiar stars with large-scale organised magnetic 
fields (Ap and Bp stars). In these stars, the abundance 
distribution of certain elements is non-uniform and non- 
symmetric with respect to the rotation axis. Numerous 
studies of Ap and Bp stars have revealed a kind of symme- 
try between the topology of the magnetic field and the ele- 
ment distribution. Thus, the structure of the magnetic field 
can be studied by measurements of the magnetic field, using 
spectral lines of each element separately. Although strong 
large-scale magnetic fields have not generally been found in 
HgMn stars, it has never been ruled out that these stars 
might have tangled magnetic fields of the order of a few 
thousand Gauss with only very weak net longitudinal com- 
ponents. Bychkov et al. (|2009j) compiled all published lon- 
gitudinal magnetic field measurements of different groups 
of chemically peculiar stars and concluded that the group 
of HgMn stars possesses the second weakest fields after the 
group of Am stars. 

In the last few years, a number of attempts to de- 
tect mean longitudinal magnetic fields in HgMn stars have 
been made by several authors using the line addition tech- 
nique, the least-squares deconvoluti on (LSD ), most recently 
by Makaganiuk et al. (e.g. I2011al 12011b! I2011c|) . A high 
level of precision, from a few to tens of Gauss, is achieved 
through application of the LSD technique (Donati et al. 
1997), which combines hundreds of spectral lines of vari- 
ous elements. This technique assumes that all spectral lines 
are identical in shape and can be described by a scaled 
mean profile. However, the lines of different elements with 
different abundance distributions across the stellar surface 
sample the magnetic field in different manners. Combining 
them as is done with the LSD technique may lead to the 
dilution of the magnetic signal or even to its (partial) can- 
cellation, if enhancements of different elements occur in re- 
gions of opposite magnetic polarities. Using this technique, 
Makaganiuk et al. (|2011b[) analysed HARPS spectra for a 
sample of HgMn and normal B-type stars and reported no 
detection at a 3<r level in any of the studied targets. 

Strangely enough, although the authors were aware of 
the inhomogeneous distribution of the elements on the sur- 
face of HgMn stars, no analysis has been done on the lines 
of individual elements separately. Only in the most recent 
Makaganiuk et al. Q2012p study did the authors decide to 
constrain the longitudinal magnetic field of HD 11753 by 
computing mean profiles for Y, Ti, and Cr, which show dif- 
ferent spot distributions on the stellar surface. Their results 
indicate that the upper limit for the strength of the mag- 
netic fields can reach about 20-30 G, with a typical error 
bar from 8 G to 15 G. 

To test if the use of the LSD technique may indeed 
account for the non-detection of magnetic fields in HgMn 
stars, we decided to re-analyse the HARPS spectropolari- 
metric material that recently became publically available 
in the ESO archive. In this work, we use a completely dif- 
ferent approach for the measurements of magnetic fields, 
namely, the moment technique developed by Mathys (e.g. 
H99T1 11995al I1995b|) . This technique allows us not only to 
determine the mean longitudinal magnetic field, but also 
to prove the presence of the crossover effect and quadratic 
magnetic fields. This information cannot be obtained from 
the LSD technique, as it assumes that all spectral lines are 
identical in shape and can be described by a scaled mean 
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Table 1. List of the studied targets. 



Object 


Other 


V 


Spectral 


Instr. 


Name 


Identifier 




Type 




HD 11753 


4>Phe 


5.11 


B9V 


HARPS 


HD 19400 


eRyi 


5.50 


B3V 


FORS2 


HD 27376 


41Eri 


3.55 


B9V, SB2 


HARPS 


HD 32964 


66 Eri 


5.10 


B9V, SB2 


HARPS 


HD 33904 


/iLep 


3.28 


B9IV 


HARPS 


HD 34364 


ARAur 


6.14 


B9V, SB2 


SOFIN 


HD 53244 


7 CMa 


4.10 


B8II 


HARPS 


HD 65949 


CPD-60 966 


8.37 


B8/B9 


FORS2 


HD 78316 


k Cnc 


5.24 


B8III, SB2 


HARPS 


HD 101189 


HR4487 


5.14 


B9IV 


HARPS 


HD 221507 


/3Scl 


4.37 


B9.5IV 


HARPS 


(Superficially) normal B-type stars 


HD 179761 


V1288Aql 


5.14 


B8II-III 


HARPS 


HD 209459 


21 Peg 


5.82 


B9.5V 


HARPS 



Notes. Spectral types and visual magnitudes are taken from 
SIMBAD. 



profile. Furthermore, we present five new spectropolarimet- 
ric observations of the eclipsing system AR Aur with a pri- 
mary star of HgMn peculiarity, obtained with the echelle 
spectrograph SOFIN at the Nordic Optical Telescope at 
the end of 2010. The variability and spot distribution of 
this star has been intensively studied by our group in the 
few last years (Hubrig et al. 2006a, 201TJ) . In addition, a 
new set of polarimetric spectra was obtained for the HgMn 
star HD 65949 and the hotter analog of HgMn stars, the 
PGa star HD 19400, using FORS2 installed at the VLT. 
All targets discussed in this work are presented in Table [1] 
together with their visual magnitudes and spectral types. 

In Sect.[5J we describe the observations, data reduction, 
and methods of our magnetic field measurements. The stel- 
lar characteristics and the magnetic field measurements for 
the individual targets are reviewed in Sect. [31 Finally, in 
Sect. SI we discuss our results in the context of possible 
mechanisms at play, which can be considered responsible 
for the development of the inhomogeneous abundance dis- 
tribution and the presence of weak magnetic fields. 

2. Observations and magnetic field measurements 

2.1. HARPS observations 

The majority of the spectra analyzed in this work were ob- 
tained with the HARPSpol polarimeter (Snik et al. [2011) 
feeding the HARPS spectrometer (Mayor et al. I2003P at 
the ESO 3.6 m telescope on La Silla. We downloaded from 
the ESO archive the publically available multi-epoch and 
single-epoch spectra for eight HgMn stars, one normal B- 
type star, and one superficially normal B-type star ob- 
tained with HARPSpol in 2010. Out of this sample, three 
HgMn stars belong to spectroscopic double-lined systems. 
The study of magnetic fields in one normal star and one 
superficially normal B-type star is of special interest be- 
cause they were reported as weakly magnetic in previous 
polarimetric studies. All spectra have a resolving power of 
R = 115 000 and a signal-to-noise ratio (S/N) between 80 
for HD 209459 and 589 for 41 Eri. The information on the 
individual observations, including the MJD dates and the 
S/N values, is given in Tabled To distinguish between SB1 
and SB2 systems, we will call the SB1 systems with their 



Table 2. Logbook of the spectropolarimetric observations. 



Object 


MJD 


S/N 


HARPS observations 


HD 11753 


55199.0975 


332 




55200.0549 


299 




55201.0296 


395 




55202.0438 


433 




55203.0356 


348 




55204.0228 


430 




55205.0518 


353 




55206.0271 


406 




55207.0207 


476 




55209.0147 


379 




55210.0368 


564 


41 Eri 


55201.2687 


507 




55210.2561 


422 




55212.2414 


589 




55213.2305 


411 


66 Eri 


55202.2262 


245 




55203.2457 


299 




55204.2384 


320 




55205.2570 


185 




55206.2375 


304 




55207.2363 


315 




55209.2426 


261 




55210.2476 


216 




55211.2523 


241 




55212.2290 


380 


HD 33904 


55204.2676 


600 


HD 53244 


55204.2799 


477 


k, Cnc 


55202.3038 


230 




55211.2873 


339 


HD 101189 


55201.3629 


179 




55201.3629 


179 


HD 179761 


55319.3267 


123 


HD 209459 


55417.0999 


80 




55417.1395 


87 




55421.2527 


257 


HD 221507 


55210.0485 


349 


SOFIN observations 


ARAur 


55544.1347 


288 




55553 2003 


310 




55554.0300 


365 




55555.1200 


342 




55556.0350 


346 


FORS 2 observations 


HD 19400 


55845.2952 


1604 




55935.1094 


1322 


HD 64949 


54108.1771 


794 




54433.3662 


876 




55686.0703 


1304 




55688.0404 


911 



HD number and use the other identifiers given in the second 
column of Table [T] for the SB2 systems. 

The HARPS archive spectra cover the wavelength range 
3780-6913 A, with a small gap around 5300 A. Each obser- 
vation of the star is usually split into four to eight sub- 
exposures, obtained with four different orientations of the 
quarter-wave retarder plate relative to the beam splitter 
of the circular polarimeter. The reduction was performed 
using the HARPS data reduction software available at the 
ESO headquarters in Germany. The Stokes / and V param- 
eters were derived following the ratio method described by 
Donati et al. (1997), ensuring in particular that all spuri- 
ous signatures are removed at first order. Null polarisation 
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Table 3. Measurements of the mean longitudinal magnetic field using sets of lines belonging to different elements. All 
quoted errors are la uncertainties. 



MJD 


Phase 


<S Z >Ti 


















[G] 


[G] [G] 


[G] 


[G] 


[G] 


[G] 


[G] 



HD 11753 



55203.0356 


0.098 


0±16 


12±16 


18±11 


-1±20 


55204.0228 


0.203 


-39±12 


0±13 -14±11 


-24±13 


-48±11 -17±13 


55205.0518 


0.308 


-3±10 


-18±13 


-7±11 


9±15 


55206.0271 


0.413 


8±9 


-11±16 


-2±11 


10±15 


55207.0207 


0.518 


1±10 


8±11 


12±10 


-3±10 


55199.0975 


0.678 


7±12 


20±15 


18±16 


4±15 


55209.0147 


0.729 


12±11 


20±8 


8±11 


18±19 


55200.0549 


0.783 


26±9 


63±19 -12±17 


25±7 10±8 


45±19 




U.OOi 




-9±16 


15±11 


44-1 1 


55201.0296 


0.888 


32±21 


14±23 


20±10 


43±20 


55202.0438 


0.993 


20±14 


22±16 


1±6 


24±15 






41 Eri A - 


- more massive, with stronger 


Hg and Fe lines 




55210.2561 


0.276 


-89±22 


31±25 


-38±23 




55201.2687 


*0.483 


-6±21 




-3±14 




55212.2414 


0.673 


36±25 




29±10 




55213.2305 


0.870 


98±47 




39±33 








41 Eri B 


- less massive with stronger Mn and Ti lines 




55210.2561 


0.276 


94±35 




81±28 




55201.2687 


*0.483 


-6±21 




-3±14 




55212.2414 


0.673 


-103±27 - 


-34±30 


-31±10 -1±11 




55213.2305 


0.870 


-72±29 




-40±22 





66 Eri A 



55203.2457 


0.063 


-17±74 




-11±43 


55209.2426 


0.194 


119±86 




32±58 


55204.2384 


0.250 


27±72 




34±26 


55210.2476 


0.384 


5±86 




-61±54 


55205.2570 


0.442 


13±74 




-1±35 


55211.2523 


*0.573 


-16±47 




-20±40 


55206.2375 


0.627 


-30±66 




-40±32 


55212.2290 


0.757 


-54±49 




-71±32 


55207.2363 


0.816 


-138±62 




-25±23 


55202.2262 


0.870 


-41±73 




-48±50 


66 Eri B 


55203.2445 


0.063 


14±38 




76±44 


55209.2415 


0.194 


-95±29 -37±31 




-34±42 


55204.2372 


0.250 


-7±32 




-13±40 


55210.2464 


0.384 


-63±40 




-6±32 


55205.2558 


0.442 


-81±63 




-41±49 


55211.2506 


*0.573 


-16±47 




-20±40 


55206.2364 


0.627 


41±35 




-30±37 


55212.2273 


0.757 


68±30 




61±33 


55207.2351 


0.816 


35±23 




24±28 


55202.2251 


0.870 


16±49 




66±43 



Notes. Our phase zero for 66 Eri refers to conjunction I. Asterisks indicate rotation phases where the spectral lines of both 
components overlap. 



spectra (labeled with n in Table [3]) were calculated by com- 
bining the sub-exposures in such a way that the polarisation 
cancels out, allowing us to verify that no spurious signals 
are present in the data. Other details of the typical acquisi- 
tion and calibration of the HARPSpol observations can be 
found in papers by Makaganiuk et al. (e.g. l20lTal 12011b)) . 

The diagnostic potential of high-resolution circularly 
polarised spectra using the moment technique has been 
discussed at length in numerous papers by Mathys (e.g. 



1993' 11995a! I1995b|) . Wavelength shifts between right- and 
left-hand side circularly polarised spectra are interpreted 
in terms of a longitudinal magnetic field (B z ). The major 
problem in the analysis of high-resolution spectra is the 
proper line identification of blend-free spectral lines. The 
quality of the selection varies strongly from star to star, 
depending on the binarity, the line broadening, and the 
richness of the spectrum. In addition, due to the inhomoge- 
neous element distribution on the stellar surface, the degree 
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MJD 


Phase 


(Bz) Ti 




(B z ) Cl (Bz) c 




(^>Fe,n 










[G] 


[G] 


[G] [G] 


[G] 


[G] 


[G] 


[G] 




HD 33904 


55204.7676 




-32±11 






-3±14 








ARAurA 


55544.1347 


0.030 


-138±97 






-160±89 




-213±108 




55553 2003 


0.223 


62±96 






-29±76 




-101±93 




55554.0300 


0.424 


266±110 






82±41 




254±88 




55555.1200 


0.687 


-199±63 


64±65 




-207±67 


33±69 


-216±68 


78±72 


55556.0350 


0.908 


-335±94 


-70±98 




-220±59 


-30±65 


-353±100 


-50±112 


ARAurB 


55544.1347 


0.030 


208±116 






263±85 


-61±80 






55553 2003 


0.223 


201±136 






140±58 








55554.0300 


0.424 


-162±61 






-116±63 








55555.1200 


0.687 


196±73 






165±74 








55556.0350 


0.908 


71±102 






93±89 








HD 53244 


55204.2799 




62±43 




85±49 


98±46 








HD 78316 


55202.3038 




-22±15 




-52±23 


-32±10 


-4±11 






55211.2873 




84±38 




-13±18 


3±17 








HD 101189 


55201.3629 




-74±24 


22±25 


-24±18 


-42±26 




-37±11 


-6±12 


HD 179761 


55319.4267 










-88±49 








HD 209459 


55417 1494 




99±22 


24±21 


68±20 -5±19 


53±17 


9±18 






55421.2598 




31±16 




12±19 


27±10 








HD 221507 


55210.0485 




78±44 




95±58 


56±24 




78±25 


-14±24 



of blending changes with the stellar rotation phase. All the 
lines we tried to employ in the diagnosis of the magnetic 
fields on the surface of our target stars are presented in 
Table SI together with their Lande factors. We note that 
the actual list of lines measured in each star on each spec- 
trum can differ from one observation to the next due to the 
different quality of the spectra and variable blending over 
the rotation cycle. For the lines of iron-peak elements, the 
Lande factors were taken from Kurucz's list of atomic data 
(Kurucz 1989), while the Lande factors for the Y lines were 
retrieved from the Vienna Atomic Line Database (VALD; 
e.g. Kupka et al.QIMI)- 

The measurements of the mean longitudinal magnetic 
field in our targets using line lists of different elements sep- 
arately are presented in Table [3] Since the element abun- 
dances are different in different stars, only the lines of ions 
with a larger number of measurable lines were used for the 
analysis. A study of elements with a small number of lines 
sets considerable limitations on the amount of information 
on the presence of a magnetic field. The most populated 
samples of lines belong to Ti and Fe, followed by Y and Cr. 
In a number of stars, Cr and Y lines appear too faint and 
could not be used in the measurements. The mean longi- 
tudinal magnetic fields were determined from null spectra 
only in the phases where the longitudinal field was detected 
at a 3er significance level. Since no significant fields could be 



determined from null spectra, we conclude that any notice- 
able spurious polarisation is absent. The measurements on 
the spectral lines of individual elements using null spectra 
are labeled with n in Table El 

To calculate the rotation phases for the star HD 11753, 
the following ephemeris was used: 



T max = 2451800.0 + (9.531 ± 0.001)E, 

where T max is the HJD for which the equivalent width of 
Yn is at maximum. 

This ephemeris is based on the work by Briquet et al. 
(|2010[) and Korhonen et al. (2012, in preparation), who 
studied numerous spectroscopic time series of this star ac- 
quired in the last decade. The discussion on the rotation pe- 
riods of the three SB2 systems, 41Eri, 66Eri, and AR Am- 
is presented in the next section, where we conclude that the 
binary components in all systems apart from the 66 Eri sys- 
tem are already synchronised and their rotation periods cor- 
respond to their orbital period. The adopted ephemeris to 
calculate rotation phases for 41 Eri and 66 Eri is discussed 
in detail in the same section, while the rotation phases for 
ARAur were calculated using ephemeris presented in the 
work of Albayrak et al. ([2"UU3")l . 
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Table 4. Line lists for the different elements used in the 
magnetic field measurements. 
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Another approach to study the presence of magnetic 
fields in upper main sequence stars is to determine the value 
of the mean quadratic magnetic field, 

{B q ) = {{B^) + {Bl)f'\ 

which is derived through the application of the moment 
technique, described, e.g. by Mathys & Hubrig (|2006[) . Here, 
(B 2 ) is the mean square magnetic field modulus, i.e. the av- 
erage over the stellar disc of the square of the modulus of 
the magnetic field vector, weighted by the local emergent 
line intensity, while (B 2 ) is the mean square longitudinal 
magnetic field, i.e. the average over the stellar disc of the 
square of the line-of-sight component of the magnetic vec- 
tor, weighted by the local emergent line intensity. 

The mean quadratic magnetic field is determined from 
the study of the second-order moments of the line profiles 



recorded in unpolarised light (that is, in the Stokes param- 

eter I). The second-order moment R} (A/) of a spectral 
line profile recorded in unpolarised light about its center of 
gravity A/ is defined as 

fif(Ao) = ^ J rjr I (A - A )(A - A ) 2 dA. 

The integration runs over the whole width of the ob- 
served line (see Mathvs 119881 for details). W\ is the line 
equivalent width; rjr { is the line profile 

t Ti = 1 - 

Tj (resp. Fi c ) is the integral over the visible stellar disk of 
the emergent intensity in the line (resp. in the neighboring 
continuum). The analysis is usually based on the consid- 
eration of samples of reasonably unblended lines and criti- 
cally depends on the number of lines that can be employed. 
Importantly, contrary to the mean longitudinal field, the 
mean quadratic magnetic field provides a measurement of 
the field strength that is fairly insensitive to its structure. 
Thus it is especially well suited to detect fields that have 
a complex structure, as is likely the case in HgMn stars. 
Using this method, Mathys & Hubrig (1995) could demon- 
strate the presence of quadratic magnetic fields in two close 
double-lined systems with HgMn primary stars, 74 Aqr and 
X Lup. The measurements of the mean quadratic magnetic 
field in our targets are presented in Table [S] 

The crossover effect can be measured by the second- 
order moment about the centre of the profiles of spectral 
lines recorded in the Stokes parameter V. Mathys (1995a) 
showed that one can derive from the measurements a quan- 
tity called the mean asymmetry of the longitudinal mag- 
netic field, which is the first moment of the component of 
the magnetic field along the line of sight, about the plane 
defined by the line of sight, and the stellar rotation axis. 
The secon- order moment of a line profile recorded in the 
Stokes parameter V with respect to the wavelength Ao of 
the centre of gravity of the corresponding unpolarised pro- 
file is defined as 

< 2) (A ) = 77^/ ^(A-A )(A-A ) 2 dA, 

where W\ is the equivalent width of the line and rjr v is 
the line profile in the Stokes parameter V. A detectable 
crossover effect was found only in very few stars. These 
measurements are discussed in the following section. 

2.2. SOFIN observations of ARAur 

Spectropolarimetric observations of the double-lined eclips- 
ing binary ARAur were obtained on five nights between 
December 12 and December 26, 2010. We used the low- 
resolution camera (R ks 30 000) of the echelle spectrograph 
SOFIN (Tuominen et al. 11999ft . mounted at the Cassegrain 
focus of the Nordic Optical Telescope (NOT). With the 
2K Loral CCD detector, we registered 40 echelle orders 
partially covering the range from 3500 to 10 000 A, with 
a length of the spectral orders of about 140 A at 5500 A. 
The polarimeter is located in front of the entrance slit of 
the spectrograph and consists of a fixed calcite beam split- 
ter aligned along the slit and a rotating super-achromatic 
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Table 5. Measurements of the quadratic magnetic field using sets of lines belonging to different elements. All quoted 
errors are 1 a uncertainties. 



MJD 


Phase 


<-Bq>Ti ( B i)cv 


<£ q > Fe 








[G] [G] 


[G] 


[G] 


HD 11753 


55206.0271 


0.413 


2034±669 






55202.0438 


0.993 


2168±502 






41 Eri A - 


more massive, with stronger Hg and Fe lines 




55213.2305 


0.870 


4520±758 






41 Eri B - 


less massive, with stronger 


Mn and Ti lines 




55212.2414 


0.673 


3885±1115 


2332±488 




55213.2305 


0.870 


4484±1201 


2106±634 




66 Eri A 


55203.2457 


0.100 




3955±1216 




55205.2570 


0.465 




4941±1174 




55206.2375 


0.642 


7687±2293 


5861±1048 




55212.2290 


0.727 


7915±1760 






55202.2262 


0.916 




5899±1940 




66 Eri B 


55203.2457 


0.100 


5480±1175 






55209.2426 


0.186 


4306±805 






55210.2476 


0.368 


4070±1239 






55205.2570 


0.465 




7271±1873 




55207.2362 


0.823 


4980±931 






55202.2262 


0.916 


3834±1247 






HD 33904 


55204.7676 




4835±1402 






ARAurA 


55553 2003 


0.223 


6339±1879 






55555.1200 


0.687 


6721±2231 






55556.0350 


0.908 


5849±1550 


8238±1750 


HD 101189 


55201.3629 




5696±709 2731±842 


3420±1041 




HD 179761 


55204.2799 






3246±1052 




HD 209459 


55417 1494 




1269±422 






55421.2598 




1470±409 







quarter-wave plate. Two spectra polarised in opposite sense 
are recorded simultaneously for each echelle order, provid- 
ing sufficient separation by the cross-dispersion prism. Two 
to four sub-exposures with the quarter-wave plate angles 
separated by 90° were used to derive circularly polarised 
spectra. A detailed description of the SOFIN spectropo- 
larimeter and its polarimetric data reduction is given in 
Ilyin (|2000p . The spectra are reduced with the 4 A soft- 
ware package (llyin 2000). Bias subtraction, master flat- 
field correction, scattered light subtraction, and weighted 
extraction of spectral orders comprise the standard steps 
of the image processing. A ThAr spectral lamp is used for 
wavelength calibration, taken before and after each target 
exposure to minimise temporal variations in the spectro- 
graph. 

The results of the longitudinal and quadratic magnetic 
field measurements using line lists for Ti and Fe are pre- 
sented in Tables |3] and [5j Due to the lower spectral reso- 
lution of our SOFIN spectra and the somewhat lower S/N 



achieved with a smaller diameter telescope, the measure- 
ment uncertainties for ARAur are on average larger than 
those achieved in the HARPS spectra. In addition, the line 
widths of the order of 23kms~ 1 in this system exceed the 
line widths for most of the targets in our sample. 

2.3. FORS1/2 measurements of the PGa star HD 19400 and 
the HgMn star HD 65949 

A few polarimetric spectra of the PGa star HD 19400 and 
the HgMn star HD 65949 were previously obtained with 
FORS 1 (Hubrig et al. I2006bl [20TT]) and most recently with 
FORS 50 on Antu (UT1) from May 2011 to January 2012. 
The star HD 19400 belongs to the group of PGa stars, pre- 
senting the hotter extension of the HgMn stars. PGa stars 
exhibit strongly overabundant P and Ga and deficient He. 



The spectropolarimetric capabilities of FORS 1 were moved 
to FORS 2 in 2009. 
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Table 6. Magnetic field measurements of HD 19400 and 
HD 65949 using FORS1/2. Already published measure- 
ments are marked by asterisks. All quoted errors are la 
uncertainties. 



MJD 


PI 


(B,) hvd [G] 


HD 19400 


* 52852.371 


151±46 


217±65 


55845.295 


14±24 


32±26 


55935.109 


-65±26 


-110±30 


HD 65949 


*53002.082 


-290±62 


-143±71 


*54108.177 


-116±32 


-131±39 


*54433.366 


-18±21 


-18±30 


55686.070 


-77±24 


-63±24 


55688.040 


-182±34 


-163±34 



Maitzen (1984) suggested the presence of a magnetic field in 
HD 19400, using observations of the A5200 feature. Similar 
to the group of HgMn stars, the spectra of HD 19400 reveal 
a strong overabundance of the elements Mn and Hg (Alonso 
et al. l2003p . This star is likely a binary system; Dommanget 
& Nys (|2002[) mention in the CCDM catalogue a nearby 
component at a separation of O'.'l and a position angle of 
179°. 

The spectrum of the HgMn star HD 65949 was studied 
in detail by Cowley et al. (120101, who discovered that this 
star exhibits enormous enhancements of the elements rhe- 
nium through mercury (Z = 75 — 80). In the catalogue of 
HgMn stars (Schneider I198ip . HD 65949 is mentioned as an 
SB1 system. 

During the observations with FORS2, we used a slit 
width of 0'.'4 and the GRISM 600B to achieve a spectral 
resolving power of about 2000. A detailed description of 
the assessment of the longitudinal magnetic field measure- 
ments using FORS 2 is presented in our previous work (e.g. 
Hubrig et al. I2004a[ I2004H and references therein). The 
mean longitudinal magnetic field, (B z ), was derived using 

V _ 9cse\ 2 Id! 

I \-nm e c? I d\ [ z) ' 

where V is the Stokes parameter that measures the circular 
polarisation, / is the intensity in the unpolarised spectrum, 
g e H is the effective Lande factor, e is the electron charge, 
A is the wavelength, m e the electron mass, c the speed of 
light, d//dA is the derivative of Stokes /, and (B z ) is the 
mean longitudinal magnetic field. 

The mean longitudinal magnetic field was measured ei- 
ther by using only the absorption hydrogen Balmer lines or 
by using the entire spectrum, including all available absorp- 
tion lines. Our measurements, acquired over a few years, 
are presented in Table [SI together with the modified Julian 
dates of mid-exposure. A few weak field detections at a sig- 
nificance level higher than 3cr were achieved for each star. 
As the rotation periods for both stars are unknown, it is 
not possible at the present stage to conclude whether the 
variability of the magnetic field is caused by a rotational 
modulation. 

With respect to the reliability of our FORS 2 mea- 
surements of rather weak longitudinal magnetic fields, we 
note that the feasibility of such measurements in stars of 



different mass and at different evolutionary stages using 
FORS 1/2 in spcctropolarimctric mode has been demon- 
strated by numerous studies during the last ten years. It 
has also been confirmed by comparing measurements of the 
magnetic field for well-studied magnetic stars with values 
from the literature. These stars have typically strong lon- 
gitudinal magnetic fields of the order of a few kG. What 
is not yet so well established are the measurements of sev- 
eral calibrators with magnetic fields well below 1000 G. The 
FORS V/I spectrum is deduced from four spectra: the two 
spectra for the ordinary and extraordinary beams coming 
out of the Wollaston prism, taken at two different settings 
of the retarder wave plate. They are combined using the 
following equation: 

v = i J \ ( r-r \ ( f°-r \ \ 

While /° 45 and /i 45 , as well as /° 45 and f+ 45 are taken 
at the same time, respectively, this is obviously not true for 
the two pairs and leads to differences in flux levels. Also, 
the transmission for the ordinary and extraordinary beams 
is not the same. In the equation used for the determina- 
tion of the V/I spectrum above, both varying transmission 
and flux levels at different times are compensated through 
the double difference, especially removing the need for a 
flatfield. Still, there might be residuals left after using the 
above equation that come from instrumental effects, e.g. 
crosstalk from linear to circular polarisation. The null spec- 
trum is usually determined as the difference between two 
consecutive V/I spectra. This computation should get rid 
of the magnetic field signatures and leave only the instru- 
mental features. Yet, different signal-to-noise ratios in the 
consecutive V/I spectra might introduce spurious null sig- 
nals. 

On the other hand, a few recent FORS 2 spectropolari- 
metric observations over the rotation/magnetic period of 
the bright classical Ap star HD 142070, with a weak mag- 
netic field, typical for a number of objects observed with 
FORS 1/2, revealed a good agreement with the measure- 
ments obtained with a high-resolution spectropolarimeter 
(Mathys et al. I2012p . These results indicate that the order 
of magnitude of the uncertainties achieved with FORS 2 is 
correct. 

3. Results for individual targets 

HD 11753: 

Among the studied targets, the star HD 11753 with 
eleven HARPS observations has the best rotation phase 
coverage. This system is an SB1 with a period of 41.489 d, 
according to the 9 th Catalogue of Spectroscopic Binary 
Orbits (Pourbaix et al. 2009). It is also an astromet- 
ric Hipparcos binary according to the CCDM catalogue 
(Dommanget & Nvs l2"rjTj2"j) . 

The spectra of HD 11753 exhibit a pronounced variabil- 
ity of Ti, Cr, Sr, and Y, and their surface inhomogeneous 
distribution was studied in detail by Briquet et al. (2010) 
and Makaganiuk et al. (|2012|) . Results of Doppler imaging 
reconstruction using data sets separated by just a couple of 
weeks revealed the evolution of chemical spots already on 
such a short time scale (Briquet et al. 2010). The Doppler 
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maps computed by Briquet et al. (2010J) and Makaganiuk 
et al. (|2012[) are roughly consistent, showing strongly over- 
abundant Y and Sr patches at the rotation phases 0.75-1.00 
and lower abundance patches around the phases 0.10-0.40. 
The Ti and Cr patches show a more complex surface distri- 
bution, although the distribution of the most overabundant 
patches roughly resembles those of Y and Sr. Furthermore, 
a significant variation in the latitudinal element distribu- 
tion and in the element abundance gradients is observed. 
Makaganiuk et al. (2012 ) assumed the time of the first night 
of observation as the zero phase, which corresponds to the 
rotation phase 0.42 in the work of Briquet et al. (2010). 
As we mention above, the rotation phases presented in 
Table [3] for the star HD 11753 were calculated following the 
ephemeris used by Briquet et al. (|2010p . 

The re-analysis of the HARPS spectra of this star 
by the moment technique using carefully selected line 
lists of the elements Ti, Cr, Fe, and Y shows remarkable 
results: The measurements on the Ti and Y lines reveal 
the presence of a weak negative magnetic field at the 
3cr significance level at rotation phase 0.203, while the 
measurements on Cr and Fe show the presence of a weak 
positive magnetic field at the 3<r significance level at 
phase 0.783. Exactly at these rotation phases become 
the patches with the lowest surface element abundances 
as well as the patches with the highest surface element 
abundances best visible on the stellar surface. Even more 
striking is the behaviour of the longitudinal magnetic 
field over the stellar surface presented in Fig. [TJ The 
measurements on all four elements indicate roughly the 
same kind of variability of the magnetic phase curve with 
predominantly negative field polarity at the location of 
the low-abundance patches and positive polarity at the 
location of the high-abundance patches. The behaviour of 
the magnetic field over the rotation cycle looks complex, 
indicating a non-sinusoidal character. It is also somewhat 
different for different elements and is very likely caused 
by differences in the surface abundance distribution for 
each element. Such a result is not unexpected since 
the lines of different elements with different abundance 
distributions across the stellar surface sample the magnetic 
field in different manners. Also, the presence of several 
clement spots with different overabundances close to each 
other as presented in Fig. 8 in the work of Makaganiuk 
et al. (|2012p for the Ti distribution indicates a more 
complex underlying magnetic field topology that cannot 
be represented by a simple dipole model. The presence 
of tangled magnetic fields on the surface of this star is 
strengthened by our 3cr detections of the mean quadratic 
magnetic field at two rotation phases, 0.413 and 0.993 
(see Table [5]) using Ti lines. As mentioned in the work of 
Makaganiuk et al. (|2012p . this element shows the most 
complex surface element distribution. In addition, a 3er 
detection for positive crossover effect, 330±108kms _1 G, 
was achieved at the phase 0.832 for the sample of Ti lines, 
and for negative crossover effect, — 467±147kms~ 1 G, at 
the phase 0.678 for the sample of Fe lines. 

41 Eri: 

With a visual magnitude of 3.55 this SB2 system is one 
of the brightest targets among our stellar sample, with 
an excellent S/N achieved in the HARPS spectra. Four 
HARPS spectra were obtained for this star, but only three 
of them can be used for measurements because one observa- 
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Fig. 1. Measurements of the mean longitudinal magnetic 
field as a function of the rotation phase for HD 11753. The 
measurements were carried out separately for the elements 
Ti, Cr, Fe, and Y (from bottom to top). Filled circles indi- 
cate 3cr measurements. 



tion was close to conjunction time, at MJD 55201.77, where 
the lines of both components overlap in the spectra. 

The atmospheric fundamental parameters for both com- 
ponents were studied by Dolk et al. (2003), who determined 
T cS = 12 750 K for the primary and T eff = 12 250 K for 
the secondary. Our inspection of the spectra belonging to 
the primary and the secondary reveal that the line profiles 
of several elements are clearly variable. Both components 
show typical HgMn peculiarities, but this is the first time 
that spectrum variability was also discovered in the spectra 
of the secondary component. 

The problem of analysing the component spectra in 
double-lined spectroscopic binaries is a tough one, but for- 
tunately several techniques for spectral disentangling have 
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Table 7. Radial velocity measurements carried out for the primary and secondary components in the SB2 system 41 Eri. 



Instr. 


MJD 


Phase 


RVa 




RV B 




FEROS 


53663. 


.2551 





.5134 


23 


.88±0. 


.13 


12. 


89±0. 


.22 


FEROS 


53664. 


.1549 





.6930 


76 


,92±0. 


.15 


-42. 


61±0. 


.25 


FEROS 


53665. 


.1027 





.8821 


60 


.21±0. 


.12 


-25. 


,57±0. 


.23 


FEROS 


53666. 


.1837 





.0979 


-18 


.26±0. 


.14 


55. 


,32±0. 


,23 


FEROS 


53668. 


.0968 





.4797 


10 


.29±0. 


.13 


26. 


,37±0. 


.22 


FEROS 


53668. 


.0983 





.4800 


10 


.27±0. 


.13 


26. 


,35±0. 


.22 


UVES 


53631. 


.2672 





.1290 


-27. 


.46±0. 


.28 


65. 


,25±0. 


.26 


UVES 


53632. 


.2467 





.3245 


-37 


.40±0. 


.25 


76. 


,06±0. 


.23 


EBASIM 


54398. 


.1707 





.1938 


-40 


.84±0. 


.24 


78. 


,41±0. 


.18 


EBASIM 


54398. 


.2361 





.2069 


-42 


.56±0. 


.24 


80. 


,13±0. 


.16 


EBASIM 


54398. 


.2990 





.2194 


-43 


.64±0. 


.26 


81. 


,42±0. 


.17 


EBASIM 


54398. 


.3462 





.2289 


-43 


.94±0. 


.29 


82. 


,05±0. 


.18 


EBASIM 


54399. 


.1231 





.3839 


-23 


.63±0. 


.26 


61. 


,37±0. 


.17 


EBASIM 


54399. 


.1944 





.3982 


-19 


.24±0. 


.24 


56. 


,87±0. 


.16 


EBASIM 


54399. 


.2501 





.4093 


-15 


.52±0. 


.24 


53. 


,16±0. 


.16 


EBASIM 


54399. 


.3325 





.4257 


-10 


.03±0 


.26 


47. 


,40±0. 


.16 


EBASIM 


54400. 


.1205 





.5830 


49 


.32±0. 


.23 


-14. 


,19±0. 


.17 


EBASIM 


54400. 


.1785 





.5946 


53 


.04±0. 


.25 


-18. 


,16±0. 


.19 


EBASIM 


54400. 


.2486 





.6086 


57. 


.43±0. 


.24 


-22. 


90±0. 


.20 


EBASIM 


54400. 


.3005 





.6189 


60 


,51±0. 


.25 


-25. 


98±0. 


.18 


EBASIM 


54401. 


.1022 





.7790 


79 


,77±0. 


.26 


-45. 


64±0. 


.17 


EBASIM 


54401. 


.2004 





.7985 


77 


.59±0. 


.27 


-43. 


67±0. 


.17 


EBASIM 


54426. 


.2828 





.8047 


77. 


12±0. 


.28 


-42. 


77±0. 


.18 


EBASIM 


54426. 


.3467 





.8174 


75. 


.34±0. 


.24 


-40. 


70±0. 


.18 


EBASIM 


54427. 


.0229 





.9524 


36 


.78±0. 


.25 


-1. 


10±0. 


.17 


EBASIM 


54427. 


.0727 





.9623 


33 


.46±0. 


.25 


3. 


,03±0. 


.16 


EBASIM 


54427. 


.1111 





.9700 


30 


.19±0. 


.23 


5. 


,83±0. 


.16 


EBASIM 


54427. 


.1729 





.9823 


25 


01±0. 


.26 


10. 


,88±0. 


.17 


EBASIM 


54427. 


.3072 





.0091 


15 


.87±0. 


.22 


20. 


,97±0. 


.17 


EBASIM 


54432. 


.2353 





.9927 


19 


.68±0. 


.21 


15. 


,90±0. 


.16 


EBASIM 


54432. 


.3133 





.0083 


15 


.72±0. 


.26 


21. 


,00±0. 


.18 


EBASIM 


54432. 


.3558 





.0168 


11 


.95±0. 


.26 


25. 


01±0. 


.17 


HARPS 


55201. 


.2686 





.4827 


12 


.25±0. 


.10 


26. 


11±0. 


.16 


HARPS 


55209. 


.0727 





.0403 


2 


,19±0. 


.12 


34. 


29±0. 


.18 


HARPS 


55210. 


.2560 





.2765 


-42 


.69±0. 


.12 


82. 


,74±0. 


.15 


HARPS 


55212. 


.2414 





.6727 


74 


,42±0. 


.12 


-38. 


33±0 


.16 


HARPS 


55213. 


.2305 





.8702 


64 


.60±0. 


.11 


-28. 


26±0. 


.16 



Notes. The HARPS observation at MJD 55209.0727 from the ESO archive had just three sub-exposures and for this reason was 
not used in the magnetic field determination. 



Table 8. Orbital and fundamental parameters for the SB2 
system 41 Eri. 



T (MJD conj i) 
Vo [kms" 1 ] 
K A [kms" 1 ] 
K B [kms" 1 ] 
a sin i [R®] 
M A sin 3 i [M©] 
M B sin 3 i [M®] 

q 

v sin ia [kms -1 ] 
v sin %b [kms^ 1 ] 
Ra sin i \R®\ 
Rb sin i [Rq] 



5.01031535±0.00001217 
54407.220102±0.0012 
18.27±0.07 
62.68±0.17 
64.70±0.15 
12.609±0.022 
0.5450±0.0029 
0.5279±0.0031 
0.9688±0.0034 



12.23±0.06 
11.78±0.07 
1.199±0.010 
1.199±0.010 



been developed in the past few years. In this work, we have 
applied the procedure of decomposition described in de- 
tail by Gonzalez & Levato (|2006|) . The resulting spectra 
were used to study the spectral variability of the compo- 
nents and the variation of radial velocities to improve the 



orbit of 41 Eri. We used all high-resolution spectra at our 
disposal obtained with FEROS, UVES, and HARPS, and 
lower-resolution (R = 20 000) EBASIM spectra. The ESO 
spectrographs UVES and FEROS are mounted on UT2 of 
the VLT at Paranal and at the 2.2 m telescope at La Silla, 
respectively. The EBASIM spectrograph is attached to the 
2.1m Jorge Sahade telescope at the CASLEO in Argentina. 
All measured radial velocities used to fit the orbit of this 
system are presented in Table [7J and the derived orbital 
parameters are listed in Table \8\ In Fig. [2j we display the 
measured radial velocities together with the calculated or- 
bital radial velocity curves. The orbit of the system 41 Eri 
is indistinguishable from a circular orbit. Thus, the solu- 
tion presented in Table [5] corresponds to an orbit fitting 
with fixed value e = 0. If e is fitted as a free parameter, we 
determine the value e = 0.0001 ± 0.0019. 

The projected equatorial rotation velocities were mea- 
sured using the mean spectrum of each component obtained 
from the disentangling of HARPS spectra. The parameter 
v sin i was calculated for several unblended Fe II and Cr n 
lines through the measurement of the position of the first 
zero of the Fourier transform of the line profiles. In these 
calculations, the appropriate linear limb darkening coeffi- 
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Fig. 2. Spectroscopic orbit of the system 41 Eri. 
Observations with FEROS are shown by red dots, with 
UVES by blue dots, with HARPS by black triangles, and 
with EBASIM by green symbols. Open symbols indicate 
the radial velocity measurements in the spectra of the sec- 
ondary component, which is the less massive star. 

cient for the corresponding stellar temperature and line 
wavelength was used. Since the orbit is circular and the 
time-scale for circularisation is longer than that for synchro- 
nisation, it can be assumed that the orbital motion and stel- 
lar rotation are already synchronised. The projected stellar 
radii and the measured v sin i-values for both components 
are presented in Table [51 

The error assigned to the stellar rotation is the standard 
error of the mean of the 25 lines measured and corresponds 
mainly to the random noise in the spectra. A more realistic 
estimate of the uncertainties should include other sources 
of error. We estimate that the influence of the instrumen- 
tal profile is small, contributing with about 0.05 km s" 1 . 
However, the limb darkening adopted for the rotational pro- 
file is not well known, since its value may differ significantly 
from the continuum limb darkening. For this reason, we as- 
sume that errors below 1% for the rotational velocities are 
unrealistic and consequently adopted a 1% error in the stel- 
lar radii as a reasonable estimate. 

In the disentangled spectra, the Hg II line at A3984 and 
the Fen lines appear to be stronger in the primary, while 
the Mn and Ti lines are stronger in the less massive star, i.e. 
in the component B. Since there is frequently confusion in 
the literature as to which component has to be considered 
as the primary and which as the secondary, we add the 
information on the spectrum appearance directly in both 
Tables [5J and [5J which present the magnetic field measure- 
ments. 

Further, we used a mass-radius (M-R) diagram to es- 
timate the inclination of the system and the absolute stel- 
lar parameters. In Fig. [3l we show the position of both 
components of 41 Eri in the M-R diagram for various pos- 
sible values of orbital inclination, together with several 
isochrones and isotherms from stellar models by Schaller 
et al. (|1992|) . Assuming T cff = 12 750 K for the primary and 
T off = 12 250 K for the secondary (see e.g. Dolk et al. [21105)1 . 
we obtain an inclination of 33.7°. The age of the system is 
estimated to be around 25-30 Myr, representing less than 
10% of its main sequence life time. 



Fig. 3. Position of the components in the system 41 Eri 
in the mass-radius diagram for inclinations in the range 
28 — 50° with a step of 2°. The red lines are isotherms 
interpolated in the stellar models. With T c g — 12 750 K 
for the primary and T c ff — 12 250 K for the secondary, we 
obtain an inclination of 33.7°. Open circles indicate the 
position of the secondary component. 

We note that the system 41 Eri actually seems to be a 
triple one: Hubrig et al. (2001) studied this system with 
diffraction-limited near-infrared ADONIS observations on 
La Silla and found a companion of K = 9.9 at a distance 
of 5'.'32 and a position angle of 162.5°. This finding was 
later confirmed by observations of Scholler et al. (2010) 
with NAOS-CONICA at the VLT. 

The variability of the spectral lines was studied in the 
high-resolution disentangled spectra for each orbital phase. 
A distinct variability is observed for Hg, Ti, and Fe lines 
in the primary and for Hg, Mn, Ga, and Fe lines in the 
secondary. The profile variations of several lines over the 
rotation cycle are presented in Fig. 2J In the upper panel, 
we present the root mean square (RMS) values of the resid- 
uals. The strongest RMS peaks correspond to the Hg line 
in the spectra of the primary and to Mn and Ga lines in 
the spectra of the secondary component. 

Since in both stars the line profiles belonging to several 
elements show significant variability, the disentangling pro- 
cess requires a dense high-resolution spectral time series. 
We note that even with the rather dense spectral time se- 
ries obtained, the results are less reliable in spectral regions 
where both components exhibit variable line profiles. The 
analysis of the Mn line-profile variations in the secondary 
indicates that Mn is concentrated in a spot located on the 
stellar hemisphere opposite to the hemisphere facing the 
companion star. Around the orbital phase 0.5, i.e. the phase 
when the stellar surface facing the primary becomes better 
visible, the Mn lines appear significantly weaker. The Ga 
lines and the Hgn line at A3984 show qualitatively a very 
similar behaviour. 

The distribution of the available measurements of the 
mean longitudinal magnetic field over the stellar surface 
in both components is presented in Fig. [5] As mentioned 
above, due to the overlap of spectral lines of components 
in one of the HARPS spectra, reliable results of the mea- 
surements can only be achieved at three rotation phases. 
3(7 detections were achieved using the sample of Ti lines in 
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Fe II A: Fe II B: Hg 




3983.2 3984.0 3984.8 4571.2 4572.0 4572.8 4628.8 4629.6 



4522.4 4523 2 3963 3 3964 ;.964 6 



4136.8 4137.6 4255.2 4256.0 



4522.4 43.25 3 



A (A) 



A (A) 



A (A) 



A (A) 



A (A) 



A (A) 



A (A) 



A (A) 



Fig. 4. Line-profile variations and the RMS values in the spectra of both components of 41 Eri over the rotation period. 
The dotted line corresponds to the mean line profiles. The scales for the normalised flux have ticks of 0.1 in the lower 
panels and 0.01 in the upper panels, respectively. 



O 



v 



v 




0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 



Fig. 5. Measurements of the mean longitudinal magnetic 
field as a function of the rotation phase for 41 Eri. The mea- 
surements were carried out separately for the elements Ti 
(bottom) and Fe (top). The solid line denotes the primary 
component, while the dashed line denotes the secondary 
component. Filled circles indicate 3cr measurements. 



the spectra of the primary at the rotation phase 0.276 as 
well as the samples of Ti and Fe lines in the spectra of the 
secondary at the phase 0.673. An interesting fact discov- 
ered in these measurements is that the stellar surfaces with 
low-abundance element spots that face the companion star 
show negative magnetic field polarity, while for the oppo- 
site hemisphere covered by high abundance element spots, 
the magnetic field is positive. 

A mean quadratic magnetic field at the 3cr level, 
4520±758G, was detected in the primary at the phase 
0.870 using the sample of Ti lines. As for the secondary, 



a quadratic magnetic field of 2332±488G and 2106±634G 
was detected at the phases 0.673 and 0.870, respectively, 
using the sample of Fe lines. Furthermore, we detect a 
quadratic magnetic field at a 3<r level, 3885±1115G at the 
phase 0.673, and a quadratic magnetic field 4484±1201G 
at the phase 0.870 using Ti lines. No crossover effect at a 
3<t level was detected for this system. 
66 Eri: 

Ten HARPS polarimetric spectra were obtained over al- 
most two orbital cycles of this SB2 system. Among them, 
one spectrum was taken close to the time of conjunction II, 
at MJD 55211.25, where the lines of both components over- 
lap in the spectra. The most recent determination of fun- 
damental and orbital parameters of 66 Eri was carried out 
by Makaganiuk et al. (|2011ap . In their work, however, there 
is a mismatch in the presentation of the orbital and funda- 
mental parameters of the system, i.e. some parameters are 
ascribed to component A, while in actuality they belong 
to component B. We re- analysed 66 Eri using all spectro- 
scopic material available to us: ten HARPS, four FEROS, 
and three UVES spectra. Furthermore, 18 spectra at a 
resolution of ~ 60 000 were obtained with the CORALIE 
echelle spectrograph, attached to the 1.2 m Leonard Euler 
telescope on La Silla in Chile. Additional radial velocity 
measurements were compiled from the literature. All these 
measurements for both the primary and the secondary com- 
ponent are presented in Tables O and [TU1 Similar to the 
study of 41 Eri, these radial velocity measurements were ob- 
tained as a result of spectral disentangling. As mentioned 
above, our phase zero for 66 Eri refers to the time of con- 
junction I, while the phase zero in the work of Makaganiuk 
et al. (|2011ap refers to time of periastron. The time of con- 
junction is more useful because it refers to the position 
with respect to the companion, while the time of perias- 
tron is generally not well defined for binary systems with 
non-zero eccentricities. Furthermore, our MJD values are 
slightly different (by about 0.0012) from those presented 
by Makaganiuk et al. (|2011b[) because they are heliocentric 
and correspond to the middle of exposures. 

In addition, we used the 83 old radial velocity measure- 
ments from Frost & Struve (|1924p . The orbital parameters 
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Table 9. Measurements of radial velocities for the primary and secondary component of the SB2 system 66Eri. 



HJD 


Phase 


Phase 


RVa 


RV B 


-2,400,000 


Conjunction 


Period 


[kms -1 ] 


[kms -1 ] 


CORALIE 


55212.5706 


0.69826 


0.97411 


127.32±0.31 


-65.82±0.23 


55212.7135 


0.72413 


0.99998 


135.30±0.31 


-73.73±0.24 


55213.5546 


0.87644 


0.15228 


111.55±0.34 


-49.47±0.23 


55213.7138 


0.90526 


0.18111 


96.41±0.32 


-34.28±0.22 


55214.5572 


0.05798 


0.33383 


6.41±0.31 


58.08±0.22 


55214.7184 


0.08717 


0.36302 


-8.67±0.30 


73.70±0.22 


55215.6575 


0.25722 


0.53306 


-59.42±0.34 


125.70±0.22 


55216.5658 


0.42169 


0.69753 


-32.08±0.34 


97.97±0.22 


55216.7108 


0.44794 


0.72379 


-20.15±0.33 


86.00±0.22 


55217.5666 


0.60291 


0.87875 


75.44±0.35 


-12.78±0.24 


55217.6999 


0.62704 


0.90289 


91.05±0.33 


-28.66±0.24 


55220.5940 


0.15109 


0.42694 


-36.04±0.35 


101.58±0.25 


55220.7149 


0.17298 


0.44883 


-43.05±0.35 


108.94±0.27 


55224.6002 


0.87651 


0.15236 


111.45±0.34 


-49.39±0.26 


55224.6936 


0.89342 


0.16927 


102.68±0.38 


-40.71±0.29 


55226.5722 


0.23359 


0.50944 


-56.75±0.33 


123.01±0.23 


55230.5862 


0.96040 


0.23627 


64.04±0.31 


-0.85±0.23 


55230.6666 


0.97496 


0.25083 


55.32±0.31 


8.38±0.23 


HARPS 


55202.7262 


0.91569 


0.19154 


91.50±0.32 


-27.64±0.18 


55203.7457 


0.10030 


0.37614 


-14.14±0.29 


80.64±0.18 


55204.7384 


0.28005 


0.55590 


-59.66±0.33 


127.41±0.19 


55205.7570 


0.46449 


0.74034 


-11.38±0.33 


78.25±0.17 


55206.7375 


0.64203 


0.91788 


100.79±0.33 


-37.08±0.18 


55207.7363 


0.82289 


0.09874 


133.31±0.33 


-70.24±0.18 


55209.7426 


0.18618 


0.46203 


-46.10±0.31 


113.46±0.18 


55210.7476 


0.36816 


0.64401 


-49.21±0.32 


116.75±0.18 


55211.7523 


0.55009 


0.82593 


40.88±0.33 


24.85±0.17 


55212.7290 


0.72694 


0.00279 


136.63±0.30 


-73.83±0.17 


UVES 


53632.8917 


0.65913 


0.93498 


109.56±1.07 


-47.27±1.29 


53662.7238 


0.06096 


0.33681 


4.35±1.22 


59.01±1.41 


53662.7254 


0.06125 


0.33710 


4.25±1.23 


59.21±1.41 


FEROS 


53663.8064 


0.25699 


0.53284 


-59.04±0.92 


126.52±1.12 


53664.6783 


0.41487 


0.69072 


-34.66±1.04 


101.19±1.15 


53665.6415 


0.58928 


0.86513 


66.72±0.91 


-3.18±1.04 


53666.7037 


0.78161 


0.05746 


139.46±0.93 


-77.86±1.13 



Table 11. Orbital and fundamental parameters for the 
SB2 system 66Eri. 



P[d] 


5.5225973±0.000002 


T (conj) 


2455208. 714±0.002 


T (per) 


2455207.191±0.010 


Vo [kms -1 ] 


32.30±0.07 


K A [kms" 1 ] 


100.36±0.18 


K B [kms" 1 ] 


102.97±0.12 


U! 


5.960±0.011 


e 


0.0833±0.0010 


a sin i [Rq] 


22.109±0.025 


M sin 3 i [Ms] 


4.761±0.016 


M A sin 3 i [M©] 


2.411±0.007 


M B sin 3 i [M©] 


2.350±0.009 


q 


0.9747±0.0022 


v sin ia [kms -1 ] 


17.68±0.04 


v sin jb [kms -1 ] 


17.00±0.05 


Ra sin i [Rq] 


1.85±0.02 


Rb sin i [Rq] 


1.78±0.02 



are shown in Table ITT1 In Fig. H2 we present our calculated 



radial velocity curves along with the old and new radial 
velocity data. 

Due to a rather large eccentricity of the system, e = 
0.0844, the rotation period of 66 Eri is not expected to be 
the same as the orbital period. The period of corotation 
at periastron is 4.657 d, while the pseudo-synchronous pe- 
riod as defined by Hut (fT9lH"j) . is 5.3015 d. The latter is 
expected to correspond to the rotation period. Makaganiuk 
et al. (|201 la[) assumed the period of the orbital motion as 
the period of spectral variability, while our rotation phases 
presented in Table [3] adopt the pseudo-synchronisation pe- 
riod. 

Similar to the procedure applied for 41 Eri, we used 
the M-R diagram to estimate the inclination of the sys- 
tem. Assuming T c g = 11 077 K for the primary and T e g = 
10 914 K for the secondary (see e.g. Makaganiuk et al. 
I2011ap . the system appears at an inclination of 73.7°. 
Figure [7] shows the position of the stellar components in 
the M-R diagram, where both companions are located very 
close to the zero age main sequence (ZAMS). 

The presence of an additional companion in this system 
was reported from diffraction-limited near-infrared obser- 
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Table 10. Measurements of radial velocities for the primary and secondary component of the SB2 system 66Eri, taken 
from the literature. 



HJD 


Phase 


Phase 


RVa 


RV B 


-2,400,000 


Conjunction 


Period 


[kms^ 1 ] 


[kms^ 1 ] 


Young (1976) 


41379.6890 


0.92001 


0.19586 


84.60±2.50 


-28.30±3.00 


41380.7080 


0.10452 


0.38037 


-15.10±2.50 


77.60±3.00 


41381.6830 


0.28107 


0.55692 


-62.30±2.50 


125.50±3.00 


41382.6700 


0.45979 


0.73564 


-17.80±2.50 


81.80±3.00 


41383.6870 


0.64394 


0.91979 


101.80±2.50 


-40.10±3.00 


42332.0070 


0.36025 


0.63610 


-52.80±2.50 


122.80±3.00 


42333.9840 


0.71823 


0.99408 


134.60±2.50 


-69.90±3.00 




Yushchenko et al. 


(1999) 




50057.3750 


0.22523 


0.50108 


-53.00±1.50 


123.10±1.50 


50060.4060 


0.77407 


0.04992 


142.90±1.50 


-74.80±1.50 




Yushchenko et al. 


(2001) 




50689.9890 


0.77532 


0.05117 


139.70±1.50 


-79.50±1.50 


50690.9820 


0.95513 


0.23098 


67.00±1.50 


-4.60±1.50 


50691.9900 


0.13765 


0.41350 


-30.50±1.50 


96.30±1.50 


50692.9860 


0.31800 


0.59385 


-59.40±1.50 


124.60±1.50 


50693.9770 


0.49745 


0.77329 


5.30±1.50 


57.70±1.50 


Catanzaro & Leto ( 200411 


52189.5966 


0.31560 


0.59145 


-64.50±3.80 


121.80±2.10 


52202.5951 


0.66929 


0.94514 


113.80±1.40 


-52.80±2.00 


52204.5894 


0.03041 


0.30626 


31.70±1.00 


31.70±1.00 


52205.5665 


0.20734 


0.48318 


-52.00±1.80 


115.40±2.20 


52220.5585 


0.92200 


0.19785 


85.40±1.80 


-25.30±1.40 


52221.5053 


0.09344 


0.36929 


-11.30±1.30 


78.00±1.50 


52235.5245 


0.63196 


0.90781 


95.20±2.60 


-34.80±0.80 


52250.4455 


0.33377 


0.60961 


-57.40±1.10 


123.10±2.20 


52566.5808 


0.57772 


0.85357 


59.60±0.90 


6.20±1.90 


52567.5663 


0.75617 


0.03202 


141.80±1.20 


-74.50±1.80 




M/M su 



Fig. 7. Position of the components of 66 Eri in the M-R di- 
agram for inclinations in the range 50—90° with a step of 2°. 
With T cS = 11077K for the primary and T cff = 10 914K 
for the secondary, we obtain an inclination of 73.7°. Both 
components are located very close to the ZAMS. Open cir- 
cles indicate the position of the secondary component. 



vations. Hubrig et al. (2001) studied this system with the 
diffraction-limited near-infrared ADONIS system on the 
ESO 3.6 m telescope and found a companion of K = 9.4 
at a distance of l'/613 and a position angle of 232.6°. The 
presence of a companion was later confirmed by observa- 
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Fig. 8. Measurements of the mean longitudinal magnetic 
field as a function of the rotation phase for 66 Eri. The 
measurements were carried out separately for the elements 
Ti and Fe. The solid line denotes the primary component, 
while the dashed line denotes the secondary component. 
Filled circles indicate 3ct measurements. 



tions of Scholler et al. (20TU|) with NAOS-CONICA at the 
VLT. 
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Fig. 6. Spectroscopic orbit of the system 66 Eri obtained using old data from the literature (left side) and new data 
consisting of ten HARPS, four FEROS, three UVES, and 18 CORALIE spectra (right side). Open symbols indicate the 
measurements of the secondary component, which is the less massive star. 



The quality of the polarimetric HARPS spectra is not 
as good as the spectra of HD 11753 and 41 Eri, with S/N 
mostly between 200 and 300. This could explain the larger 
inaccuracies in the magnetic field determinations and the 
fact that in our study only one detection at a 3ct signifi- 
cance level was achieved at the rotation phase 0.186 using 
Ti lines. The distribution of the mean longitudinal magnetic 
field values measured in both components over the rota- 
tion cycle is presented in Fig. [5J A comparison of the field 
distribution with the distribution of elements over the stel- 
lar surface presented in Fig. 7 in the work of Makaganiuk 
et al. (|2011a|) confirms the pattern already discovered in 
HD 11753 and 41 Eri: A negative mean longitudinal mag- 
netic field is measured at the location of the lower abun- 
dance patches, i.e. on the stellar surface facing the com- 
panion, and the positive mean longitudinal magnetic field 
roughly corresponds to the location of the high-abundance 
patches located on the opposite hemisphere. A positive 
3cr crossover, 1116±374kms _1 G, is observed in the sam- 
ple of Fe lines in the primary at the phase 0.384. Using 
Fe lines in the secondary, we obtain a positive crossover 
1795±396kms _1 G, at the phase 0.442. Furthermore, nega- 
tive crossover, — 1382±450kms~ 1 G, is observed in Fe lines 
in the primary at the phase 0.816. Interestingly, no sig- 
nificant crossover effect was found either in the primary 
or in the secondary using the sample of Ti lines. A mean 
quadratic magnetic field at the 3cr significance level was 
discovered in both components at several rotation phases 
using samples of Ti and Fe lines. 

HD 33904: 

Until now, four polarimetric observations have been ob- 
tained with HARPS for this HgMn star, but only one obser- 
vation was publically available in the ESO archive at the 
time of our visitor stay at the ESO headquarters, which 
was aimed at the reduction of the HARPS spectra with 



the HARPS pipeline machine. The spectral variability of 
several elements was studied by Kochukhov et al. p011[) . 
who also used the LSD technique with the mask covering 
526 lines to determine the longitudinal magnetic field and 
concluded that the upper limit for the field in this star is 
only about 3G. The available radial velocity data are not 
sufficient to decide whether this star is a member of a SB 
system, but recent adaptive optics observations by Scholler 
et al. (|2010[) presented the first direct detection of a close 
companion candidate at a separation of 0'.'352 and a posi- 
tion angle of 250.9°. Measurements using the sample of Ti 
lines reveal that the mean quadratic magnetic field is def- 
initely present on the surface of this star, while the mean 
longitudinal magnetic field measured on the same lines ac- 
counts for a 2.9<7 detection. 
ARAur: 

The eclipsing system ARAur (HD 34364, B9V+B9.5V) 
with an orbital period of 4.13 d at an age of only 4xl0 6 yr 
belongs to the Aur OBI association. Since its primary 
star of HgMn peculiarity is exactly on the ZAMS while 
the secondary is still contracting towards the ZAMS (e.g. 
Nordstrom & Johansen ll994p . it presents the best case to 
study evolutionary aspects of the chemical peculiarity phe- 
nomenon. A presence of a third body in the system was 
discovered by Chochol et al. (|1988p . The existence of the as 
yet unseen third star with a mass of at least 0.51 M was in- 
ferred from a light-time effect in the observed minima with 
a period of 23.7 yr (Albayrak et al. 120031 see also Mikulasek 
et al. I2010[) . We re-determined the spectroscopic orbit of 
this triple system using all spectroscopic data available to 
us (21 SOFIN, 20 TLS, 26 STELLA, and nine UVES spec- 
tra) . We also used the radial velocities measured by Folsom 
et al. (f20T0j) for three ESPaDOnS and four NARVAL spec- 
tra. In these calculations, we adopted the binary period and 
parameters of the light-time orbit published by Albayrak et 
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Table 12. Orbital and fundamental parameters for the 
eclipsing system ARAur. 
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Fig. 9. Spectroscopic orbit of the eclipsing system ARAur 
using SOFIN spectra obtained in 2002 and 2010 (blue 
squares), UVES spectra from 2005 (black dots), TLS- 
STELLA spectra from 2008/2009 (red triangles and green 
diamonds), three ESPaDOnS spectra (pluses), and four 
NARVAL spectra (stars). Open symbols indicate the mea- 
surements of the secondary components, which is the less 
massive star. The upper panel shows the velocity of the 
components of the eclipsing pair with respect to the center- 
of-mass of the binary; the lower panel shows the systemic 
velocity of the pair AB with respect to the center-of-mass 
of the triple system. 



al. ([2003) in order to account for the long-term variation 
of the center-of-mass velocity of the eclipsing pair and the 
light-time effect corrections to the epochs of observations. 
Radial velocity curves indicate that the eccentricity is indis- 
tinguishable from zero (e = 0.0009 ± 0.0011), in agreement 
with the photometric analysis of Nordstrom & Johansen 
(1994). We therefore adopted a circular orbit in our final 
calculations. The orbital parameters are listed in Table [T2l 
and the radial velocity curve is shown in Fig. It is clear 
from the lower panel of Fig. [9]that the influence of the third 
body cannot be neglected in the orbit fitting. 

The middle block of Table H2l corresponds to the results 
of the measurements of projected rotational velocities and 
the lower block to the combination of our spectroscopic 
orbit with the photometric inclination and relative radii 
of Nordstrom & Johansen (j 19941). According to these au- 
thors, the primary star, which is more massive and hotter, 
is slightly smaller than the secondary. However, the spec- 



Fid] 
T (conj) 
Vo [kms" 1 ] 
K A [kms" 1 ] 
K B [kms" 1 ] 
a sin i [Rq] 
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Notes. System velocity V Q corresponds to the center-of-mass 
velocity of the triple system. 



tral lines of the secondary appear slightly narrower, which, 
assuming corotation, would suggest a smaller radius. The 
difference in wsinj between both stars, however, is only 
marginal and, if real, could be explained by a non-uniform 
distribution of the elements Fe, Cr, and Ti with latitude. 
Variability of spectral lines associated with a large number 
of chemical elements was reported for the first time for the 
primary component of this eclipsing binary by Hubrig et al. 
(2006a|). In the more recent study of this system by Hubrig 
et al. (|2010p . we presented the results of Doppler imag- 
ing for the reconstruction of the distribution of Fe and Y 
over the surface of the primary. We used a spectroscopic 
time series obtained in 2005 and from 2008 October to 
2009 February. In the disentangling process, we adopted the 
photometric flux ratio obtained by Nordstrom & Johansen 
(|1994p . i.e. the light contributions of stars A and B to the 
continuum are adopted as 0.526 and 0.474, respectively. 
Our results showed a remarkable evolution of the elemental 
spot distribution and the overabundances. Measurements of 
the magnetic field with the moment technique using several 
elements revealed the presence of a longitudinal magnetic 
field of the order of a few hundred Gauss in both stellar 
components as well as a quadratic field of the order of 8 kG 
on the surface of the primary star. 

To study the magnetic field geometry on the surface 
of the components in this system and its correlation with 
the surface inhomogeneous element distribution, a number 
of nights was allocated at the NOT in December 2010. 
However, due to bad weather conditions only five polari- 
metric observations were acquired for ARAur. The distri- 
bution of the observations was fortunately rather random 
over the rotation/orbital period, allowing us to get an idea 
about the surface element distribution. The orbit of the sys- 
tem using all spectra of AR Aur available to us is presented 
in Fig. [9] The rotation phases were calculated using the 
ephemeris taken from the study of Albayrak et al. ([2003). 

In our analysis, similar to the analysis presented by 
Hubrig et al. (|2010|) . we used the improved Doppler imag- 
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Fig. 10. Maps of the abundance distribution for Fe (left), Sr (middle), and Y (right) on the surface of the primary in 
the system AR Aur. 



ing code introduced by Freyhammer et al. (|2009p . This code 
utilises Tikhonov regularisation with a grid of 6 x 6° in a 
way similar to the Doppler imaging method described by 
Piskunov (|2008p . In the reconstruction, we searched for the 
minimum of the regularised discrepancy function, which in- 
cludes the regularisation function and the discrepancy func- 
tion describing the difference between observed and cal- 
culated line profiles. All atomic data in our analysis were 
taken from the VALD data base (Kupka et al. [1999 ) . For 
the DI reconstruction, we selected the following lines be- 
longing to the elements Fe, Sr, and Y: Fen 4923.9 A, Srii 
4215.5 A, and Yn 4900.1 A. The computed Doppler maps 
in Mercator and spherical projections based only on five 
different rotation phases are presented in Fig. [TU] 

Clearly, due to the low number of observations these 
maps cannot be used to study the evolution of element dis- 
tribution over the stellar surface, but they allow us a rough 
estimate of the location of regions with element overabun- 
dances. Similar to previous maps obtained for ARAur, we 
again find that the regions with lower abundances of Y and 
Sr are located on the stellar surface facing the compan- 
ion, while spots with Y and Sr high abundances appear 
on the opposite hemispheres. Presently we have at our dis- 
posal four spectroscopic data sets obtained between 2002 
and December 2010. The temporal evolution of chemical in- 
homogeneities is presented in Fig. 111! using as an example 
the profile variations of the Srii 4215.5 A line. All four sets 
contain spectra at almost the same rotation phases, 0.68. 
The line profiles of Srii 4215.5 A observed at this phase 
in all sets show that the Sr spots definitely changed their 
shape and abundance with time (see Fig. [TTj) . The lowest 
Sr abundance appeared to be in December 2010, while the 
strongest Sr overabundance was observed in 2005. 

In our magnetic field analysis, the longitudinal magnetic 
field strength measured in AR Aur is larger than that mea- 
sured in other targets. However, due to the much lower res- 
olution of the SOFIN spectra compared with the HARPS 
spectra and the not very high S /N, our measurement accu- 
racies are poorest among the values presented in Tabled At 
present, it is not clear to us whether the stronger magnetic 
field in this system can be explained by the aspect effect, 
if we assume that the field topology is characterised by the 
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Fig. 11. Line profiles of the Srii 4215.5 A line around ro- 
tation phase 0.68 for all available data sets between 2002 
and 2010. The thick red continuous line refers to observa- 
tions with SOFIN in 2010, the blue continuous line to those 
with UVES in 2005, the dashed line corresponds to observa- 
tions at the Thiiringer Landessternwarte and STELLA, and 
the dotted line corresponds to observations with SOFIN in 
2002. 



presence of magnetic regions predominantly in the vicinity 
of stellar equators. Since the ARAur system is an eclips- 
ing binary, it is the only system in our sample that offers 
the observer the best visibility of the surface equatorial re- 
gions. In any case, this system is one of the most promising 
targets to constrain the magnetic field topology in bina- 
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Fig. 12. Measurements of the mean longitudinal mag- 
netic field presented as a function of the rotation phase 
for ARAur. They were carried out separately for the ele- 
ments Ti, Fe, and Y (from bottom to top). The solid line 
denotes the primary component, while the dashed line de- 
notes the secondary component. Filled circles indicate 3cr 
measurements. 



ries with late B-type primaries. 3<r longitudinal magnetic 
field detections were achieved for measurements using Ti, 
Fe, and Y lines in the spectra of the primary in the phases 
0.687 and 0.908 as well as for measurements using Fe lines 
in the spectra of the secondary at the phase 0.030. The 
variations of the mean longitudinal magnetic field in both 
components over the rotation cycle presented in Fig. 
display a characteristic behaviour similar to that found for 
the binary stars discussed above. The mean quadratic mag- 
netic field is detected at three different rotation phases in 
the primary. 
HD 53244: 

The SB nature of this target was suggested by Schneider 
(1981). Scholler et al. (|2010p found a companion candidate 
to this star at a separation of 0'.'332 and a position angle 
of 114.8°. The rotation perio d of H P 53244, P = 6.16d, 
was derived by Briquet et al. (2010) from the study of the 
variations of radial velocities and equivalent widths of spec- 
tral lines belonging to the inhomogeneously distributed ele- 
ments Hg and Mn. However, Y lines appear extremely weak. 
No mean longitudinal magnetic field was detected using the 
single HARPS observation of this star. A negative crossover 
effect at the 3a significance level, — 4073±1118kms~ 1 G, 
was detected using the sample of Cr lines. 



k Cnc: 

This system is an SB2 with a period of 6.3933 d, accord- 
ing to the 9 th Catalogue of Spectroscopic Binary Orbits 
(Pourbaix et al. EffiM)) . Mason et al. pDUTj) found a com- 
panion on Besselian year 1999.1606 at a separation of 0'.'286 
and a position angle of 108°. This star was also resolved by 
Roberts et al. (|2005|) on Besselian year 2003.0077 with a 
separation of 0'.'30 and a position angle of 104°. Scholler et 
al. ( 2010) detected this companion at a separation of 0'.'269 
and a position angle of 109.7°. 

Two HARPS polarimetric spectra of this system are 
available in the ESO archive. The detection of very faint 
spectral lines of the secondary in both spectra allowed us 
to determine the mass ratio, q = 0.490 ± 0.003. Combining 
our spectroscopic mass ratio with the orbital parameters P, 
e, and Ka of the published SB1 orbit by Aikman ([1976), 
we obtain M A sin 3 i = 3.74 ± 0.13 M Q and M B sin 3 i = 
1.83 ± 0.04 Mq. To learn about the nature of the secondary 
component and measure its rotational velocity, we applied 
spectral disentangling. The obtained secondary spectrum 
is consistent with a mid-A spectral type, which contributes 
only 7% to the total light of the system at A5000. Since 
most lines in the spectrum of the secondary component 
are blended, the rotation velocity was determined using the 
technique developed by Diaz et al. (2011). As a result, we 
obtain vsini = 42.3±1.5kms~ 1 . The primary star presents 
very narrow spectral lines with the v sin i- value below 5 
kms -1 . Both components rotate asynchronously. In fact, 
according to estimates of their radii, the usini values for 
pseudo-synchronisation for stars A and B are about 33 and 
14 km s , respectively, and thus very far from the observed 
values. 

With an orbital period of 6.39 d, the orbital phase dif- 
ference between the two HARPS observations is 0.41. The 
strongest lines in the spectra belong to the elements Mn and 
Fe, while Y lines are extremely weak. The weak negative 
longitudinal magnetic field is detected at the 3cr significance 
level in the HARPS spectrum obtained at MJD 55202.304 
using Fe lines. Also, measurements using Mn lines result- 
ing in (B z ) = —47 ± 15 G confirm the presence of a weak 
negative magnetic field at this epoch. No significant detec- 
tion was achieved for the crossover effect and the quadratic 
magnetic field. 

HD 101189: 

This target is not known to belong to a SB system. 
Scholler et al. (2010) found a close companion candidate to 
HD 101189 at a separation of 0'.'337 and a position angle of 
104.1°. 

The first HARPS spectropolarimetric observation of 
this star was carried out at MJD 54982.4744, and its anal- 
ysis was reported by Hubrig et al. (|201ip . For the same 
HARPS observations, Makaganiuk et al. (2011b) reported a 
non-detection with (B z ) = —95 ± 66 G using the LSD tech- 
nique. The quality of the spectra obtained at this epoch was 
especially poor, with achieved S/N between 12 and 52. The 
measurements of Hubrig et al. delivered the same magnetic 
field (of the order of — 200 G) for both the regular science 
Stokes V observations and the null spectrum. As the au- 
thors showed in their Fig. 8, this is because the spectrum 
with the best S/N ratio (52 in that case) significantly con- 
tributes to the null spectrum. This result implies that if 
the observations are carried out with low S/N, and the dif- 
ference in S/N in the sub-exposures is large, the magnetic 
field cannot be measured conclusively, and there is no use 
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in the null spectrum to prove whether or not the detected 
field is real. 

The second HARPS observation of HD 101189 was ob- 
tained on MJD 55201.3629. For this observation, we detect 
a weak magnetic field at the 3a level (B z ) = — 74 ± 24 G 
using Ti lines and (B z ) = —37 ±11 G using Y lines. A 
quadratic field at 3er significance level was detected for all 
line samples apart from the sample of Y lines. The shape 
of the line profiles of different elements is very different: Y 
lines appear double with two absorption maxima in the line 
profiles, while this profile shape is not observed in Mn, Fe, 
Cr, and Ti lines, suggesting a different spot distribution of 
these elements over the stellar surface. Only for the sample 
of Y lines do we observe the presence of a weak positive 
crossover effect, 764±338kms" 1 G, at the 2.3ct level. 

HD 221507: This target is not known to belong to a 
SB system. Hubrig et al. (2001) studied this target with 
ADONIS and did not detect any companion. However, 
Scholler et al. (|2010p found a companion candidate to 
HD 221507 at a separation of 0'.'641 and a position angle 
of 240.2°. 

The rotation period of this star, P = 1.93 d was derived 
by Briquet et al. (2010) from the study of the variations 
of radial velocities and equivalent widths of spectral lines 
belonging to the inhomogeneously distributed elements Hg, 
Mn, and Y. We detect a weak longitudinal magnetic field 
at 3a level {B z ) = 78 ± 25 G using Y lines. Also, measure- 
ments using Mn lines, (B z ) = 75±23G, show the pres- 
ence of a positive longitudinal magnetic field. No significant 
crossover effect and quadratic magnetic field are detected 
in our measurements. 

HD 179761 and HD 209459: 

The star HD 179761 is considered in the literature as a 
normal B-type star and HD 208459 as a superficially normal 
B-type star. Superficially normal B-type stars are generally 
indistinguishable from the normal stars in terms of their 
iron-peak elemental abundances (e.g. Smith & Dworetsky 
1993). They are usually sharp- lined stars that possess nor- 
mal MK classifications, but nonetheless exhibit mild pecu- 
liarities when observed at high spectral resolution. 

Although HD 209459 is not a typical HgMn star, it is 
listed in the catalogue of HgMn stars by Schneider (|198ip . 
Our analysis of the HARPS spectropolarimetric material 
for this star indicates the absence of Hg and Y lines and 
the presence of weak Mn lines. The lines of Ti, Cr, and Fe 
appear rather strong and reveal a weak variability. 

Ti, Cr, Mn, and Y lines are weak in the HARPS spec- 
trum of HD 179761. We cannot report anything on the vari- 
ability of this star as only a single HARPS observation is 
available. 

Hubrig et al. (|1999b|) and Hubrig & Castelli ([MJT]) used 
an anomalous strength of the Fen A6147.7 line relative to 
the Fen A6149.2 line in stars with magnetic fields to study 
the presence of magnetic fields in HgMn, normal, and super- 
ficially normal late-B type stars. In this method, previously 
introduced by Mathys (|1990[) . the observed relative differ- 
ences between the equivalent widths of the two Fen lines 
are compared with those derived from synthetic spectra 
computed by neglecting magnetic field effects. While the 
differences between the equivalent widths for HD 209459 
were found to lie within the error limits, the differences for 
HD 179761 were found to be much larger, even if all ob- 
servational and computational uncertainties are taken into 
account (Hubrig & Castelli 2001), which suggests the pres- 



ence of a magnetic field in this star. Both HD 179761 and 
HD 209459 we re also observed with FORS 1 at the VLT by 
Hubrig et al. (I2006b| , who reported a 3cr detection for the 
mean longitudinal magnetic field of HD 179761. 

The observations with HARPS indicate the presence 
of a weak longitudinal magnetic field in HD 209459 at 
MJD=55417.149, measured using the samples of Ti, Cr, 
and Fe lines. A quadratic magnetic field of the order of 
3kG at the 3a level was detected in HD 179761 using Fe 
lines, and much weaker quadratic fields of the order of 1.3- 
1.5 kG were detected in both observations of HD 209459 
using Cr lines. The presence of a ^3 kG quadratic field in 
HD 179761 agrees well with the previous analysis of the 
relative magnetic intensification of the two Fell lines by 
Hubrig & Castelli (|2001|) . We note that it is not the first 
time that a quadratic magnetic field is detected in upper 
main sequence stars considered as normal or superficially 
normal stars. Mathys & Hubrig (2006) studied the presence 
of a quadratic magnetic field in the superficially normal A- 
type star HD 91375. Their analysis led to the detection of 
a 1.7 kG quadratic field at a significance level of more than 
4a. 

4. Discussion 

The region of the main sequence centred on A and B stars, 
also referred to as the "tepid stars", represents an ideal 
laboratory to study a wide variety of physical processes 
that are at work to a greater or lesser extent in most 
stellar types. These processes include radiation driven dif- 
fusion, differential gravitational settling, grain accretion, 
magnetic fields and non-radial pulsations. Understanding 
them is becoming increasingly important for the refine- 
ment of stellar evolution models and for the improved treat- 
ment of the stellar contribution in studies of galactic evolu- 
tion. Presently, the radiative diffusion hypothesis developed 
by Michaud (|1970p is the most frequently accepted the- 
ory for producing the observed anomalies in HgMn stars. 
According to Michaud et al. (|1974) . the effect of diffusion 
would cause a cloud of mercury to form high up in the 
HgMn star atmospheres. However, the impact of the ra- 
diative diffusion process has not yet been studied for any 
element considered in our work. 

Based mostly on ESO/HARPS archive data, our analy- 
sis suggests the existence of intriguing correlations between 
magnetic field, abundance anomalies, and binary proper- 
ties. In the SB2 systems with the synchronously rotating 
components, 41Eri and ARAur, the stellar surfaces fac- 
ing the companion star usually display low-abundance ele- 
ment spots and negative magnetic field polarity. The sur- 
face of the opposite hemisphere, as a rule, is covered by 
high-abundance element spots and the magnetic field is pos- 
itive at the rotation phases of the best-spot visibility. Also, 
Measurement results for the SB1 system HD 11753 indi- 
cate that element underabundance (respectively overabun- 
dance) is observed where the polarity of the magnetic field 
is negative (respectively positive). In the case of the very 
young pseudosynchronously rotating system 66Eri, addi- 
tional future observations are urgently needed to determine 
the period of the spectral variability and to investigate the 
presence of element-spot evolution. It is quite possible that 
because the system is not yet synchronised, a dynamical 
spot evolution is more pronounced compared to synchro- 
nised systems. 
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Fig. 13. Circularisation times for 41Eri (dotted lines), 
66Eri (solid lines), and ARAur (dashed lines), assuming 
that they formed with eccentric orbits with eo = 0.2 and 
0.5. 



Although only rather weak longitudinal magnetic fields 
are detected on the surface of our targets, the numerous 3er 
detections of mean quadratic magnetic fields strongly sug- 
gest that magnetic fields are present in their atmospheres. 
We note again the advantage of carrying out quadratic mag- 
netic field measurements as they provide measurements of 
the magnetic field strength that are somewhat insensitive 
to the magnetic field structure, and especially suited for 
magnetic fields with complex structure. 

The interactions of differential rotation, magnetic fields, 
and meridional flows in a tidal potential may lead to 
processes that can explain complex surface patterns. 
Possible instabilities should be studied by non-linear, three- 
dimensional simulations. The resulting fluctuations in ve- 
locity and magnetic field may not only explain the relatively 
small scales of surface magnetic fields, but also provide en- 
hanced transport of angular momentum that explain the 
high rates of synchronisation. Synchronisation in an orbit 
with about 5 d period, an initial rotation period of 1 d , and 
stellar parameters similar to 66 Eri is roughly 10 4 times the 
dissipation time scale in the star fZahn [2"008p . In a non- 
turbulent, non-magnetic star where only the microscopic 
gas viscosity can couple the surface with the interior, it 
would be impossible to reach synchronisation within the 
age of the universe. There are also gravity waves that lead 
to synchronisation, but for stars of 2-3 solar masses, this 
effect still takes roughly 10 9 years. Large-scale magnetic 
fields or turbulence can provide enough coupling to reach 
(pseudo-)synchronisation in about a million years. 

As an illustration, we calculated the circularisation 
times for the systems 41 Eri, 66 Eri, and ARAur, assuming 
that they formed with eccentric orbits. With the present 
parameters of the stars, the circularisation time scales, 
(— dlne/dt) -1 could be as high as 10 10 yr. As they evolve 
within the main sequence, circularisation becomes more ef- 
ficient, but even so they will arrive at the terminal age main 
sequence (TAMS) with a significant fraction of their origi- 
nal eccentricity: e/eo=0.15-0.70, depending on the original 
eccentricity eo. Figure [13] shows the evolution of e for the 
three systems, starting from eo = 0.2 and 0.5. The solid 
line corresponds to 66 Eri, the dotted line to 41 Eri, and 



the dashed line to ARAur. The vertical end of the curves 
to the right mark the TAMS. The calculations were per- 
formed using the stellar models calculated by Claret (2004), 
which include the internal structure constants required to 
calculate the orbital evolution. Considering that all three 
systems are young, with ages well below 10 8 yr, it would be 
expected that their present eccentricities woould not differ 
significantly from the original ones. The observed very low 
or null eccentricities cannot be explained as a consequence 
of standard tidal evolution. 

The magneto-rotational instability is a candidate that 
can provide turbulence for a wide range of magnetic field 
strengths and field topologies, but only if the internal an- 
gular velocity of the star decreases with the distance from 
the rotation axis. Since synchronisation in the systems con- 
sidered here is in fact a braking process, the star may ex- 
hibit a slower rotation near its surface than deeply inside. 
Numerical simulations should elucidate the expected en- 
hancement of the viscosity in the star and will also deliver 
surface topologies of the fields resulting from this or other 
instabilities that can be compared to the observed patterns. 

Tidal forces strongly depend on the distance between 
the components. For the three SB2 systems with synchro- 
nised and pseudosynchronised components, 41 Eri, 66 Eri, 
and ARAur, we estimated the fractional radii rj — Rj/a. 
Assuming for 41 Eri 



(v sin i) A 
{v sin i) b 
K A + K B 

e 



12.23 ± 0.06 km/s, 
11.78 ± 0.07 km/s, 
127.38 ±0.23 km/s, 
0, 



it follows that R A /r = 0.096 ± 0.001 and R B /r = 0.093 ± 
0.001 

If we assume for 66 Eri pseudosynchronisation (P ro t = 
5.3015 d), we obtain R A /r = 0.084 ± 0.001 and R B /r = 
0.080 ±0.001. 

For ARAur we use 



(v sin i) A 
(v sin i) b 
K A + K B 



22.87 ±0.08 km/s, 
22.35 ±0.07km/s, 
225.28 ±0.18 km/s, 




and obtain R A /a = 0.102±0.001 and R B /a = 0.099±0.001. 

Since the fractional radii in these systems are quite 
large, the components in the close binary systems are sub- 
ject to tidal forces acting differentially throughout their 
bodies and causing various effects in the system. The tidal 
effects include synchronisation of the rotational periods 
with the orbital period, circularisation of the orbit, and 
alignment of the rotation axes with the orbital axis. All 
these changes in the dynamical properties of the stars hap- 
pen on different time scales. The spin-down may be a sur- 
face effect, but the star then possesses an internal differ- 
ential rotation that may lead to hydrodynamic or magne- 
tohydrodynamic instabilities. Even a moderately turbulent 
state generated by an instability will lead to much stronger 
coupling with the stellar interior and redistribute angular 
momentum. 
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Fig. 14. H-R diagram for the HgMn stars in our sample. Continuous lines present the ZAMS and TAMS. Dotted lines 
present isochrones for log(t) from 7.00 to 8.70. 



To establish the evolutionary status of the HgMn stars 
in our sample, we studied their distribution in the H-R di- 
agram. To this aim, the distances were calculated mainly 
using Hipparcos parallaxes (van Leeuwen I2007p . For stars 
with high galactic latitudes, we estimated the interstellar 
absorption using Schlegel et al. (|1998[) maps and the dis- 
tance by applying the same procedure as Bilir et al. (2008). 
Since this approach is not valid for targets located very 
close to the galactic plane, we estimated the absorption in 
a different way for the low-latitude stars HD 101189 and 
ARAur.. HD 101189 is located in the field of the open 
cluster NGC3766, but it is a foreground star. This fact 
allowed us to estimate the absorption to the star using 
cluster reddening and the distance ratio. For the eclips- 
ing binary ARAur, we adopted the absolute stellar pa- 
rameters derived by Nordstrom & Johansen (|1994p from 
light and radial velocity curves. The star HD 65949, which 
is the faintest and most distant star, belongs to the open 
cluster NGC2516. We used, therefore, the cluster distance 
and reddening published by Terndrup et al. (2002) to ob- 
tain absolute magnitude and intrinsic color. In the case of 
double-lined spectroscopic binaries, the position of the in- 
dividual components in the H-R diagram was calculated 
using the spectroscopic mass ratio, from which the magni- 
tude difference between the components was estimated by 
interpolating in the Geneva stellar models (Schaller et al. 
1992). FiaurefMl shows the position of the 11 HgMn stars 
in the color-magnitude diagram, together with isochrones 
taken from Schaller et al. (|1992[) . Six stars in our sam- 
ple, 41EriA, 41EriB, 66EriA, ARAur A, HD 65949, and 
HD 221507, are located very close to the ZAMS. 

From the knowledge of orbital parameters, fundamental 
parameters of components, and fractional radii, it is possi- 



ble to estimate the reflection effect, gravitational distortion, 
and darkening for a typical HgMn star in a SB2 system. 
One of the consequences of the proximity of the stars is 
the departure from sphericity of the stellar surface. The 
star acquires a somewhat larger dimension along the line 
passing through the two star centres, while it has a smaller 
radius in the direction that is perpendicular to the former 
line and also lies in the orbital plane. The relative differ- 
ences between both radii are 0.09% and 0.08% for the two 
components of the system 41 Eri, 0.05% and 0.06% for the 
system 66 Eri, and 0.09% and 0.10% for the system ARAur. 
Observational effects of such stellar shapes are, for exam- 
ple, the variation of the ushu-value, which in the three 
systems studied is below the observational errors, and light 
variations due to the variation of the stellar disk size, which 
is expected to be of the order of 2-3mmag. An effect that 
could affect physical conditions in the stellar atmosphere 
is the reflection effect. Assuming an albedo equal to 1 (the 
usual assumption for stars with radiative envelopes) and the 
temperatures and geometrical configuration of each system, 
we estimate that the temperature excess in the surface re- 
gion facing the companion is between 20 and 40 K, which 
represents 0.18-0.38% of the star temperature. The ques- 
tion to solve is whether such small temperature differences 
can be considered a reasonable explanation of the origin of 
the low- and high-abundance spots observed either on the 
surface facing the companion or on the surface hidden to 
the companion. Although it is currently not easy to give an 
answer without a specific model of time-dependent diffu- 
sion processes, it seems very unlikely that these processes 
are the main cause of the observed abundance pattern. We 
note that the early circularisation and a very specific con- 
figuration of chemical spots in binary systems support the 
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idea of an important impact of magnetic fields on the phys- 
ical processes taking place in these stars. Since the orbital 
parameters of close binaries can change as a consequence 
of the action of tidal forces or other mechanisms, the influ- 
ence of the presence of a stellar companion on the physical 
processes playing a role in the formation of chemical anoma- 
lies is expected to vary with time. The atmospheric physical 
conditions for the development of chemical inhomogeneities 
might appear favourable at certain stages of the stellar life 
and unfavourable at other stages. 

Further studies using high S/N circularly polarised spec- 
tra well-distributed throughout the rotation period of a 
representative sample of HgMn spectroscopic binaries sam- 
pling a wide range of inclinations of the rotation axis to the 
line of sight are urgently needed to map their magnetic field 
and elemental abundances with a Zeeman-Doppler imaging 
code. The sampling of a wide range of inclination of the 
rotation axis to the line of sight in HgMn systems is es- 
pecially important because the latitudinal information on 
the surface abundance distribution of different elements in 
HgMn stars is still rather poor. In this way, it will be pos- 
sible to probe correlations between the binary properties, 
the magnetic field structure, and the abundance inhomo- 
geneities. At the same time, we will be able to gain insight 
into the mechanisms responsible for the chemical inhomo- 
geneities and their dynamical evolution on the surface of 
HgMn stars. 
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